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I. INTRODUCTION 

It is generally agreed that gravitational energy exists, but because of the equivalence principle it cannot be localized 
[1], A suitable expression which could provide a physical reasonable description of the energy-momentum density 
for gravitating system has lone been sought. All candidates had several shortcoming. In particular they violated a 
fundamental theoretical requirement—that gravitational energy should be positive—as well as requirements concerning 
localization and reference frame independence. 

During the last 30 years one of the greatest achievements in classical general relativity is certainly the proof of 
the positivity of the total gravitational energy, both at spatial and null infinity [2-4]. It is precisely positivity that 
makes this notion not only important (because of its theoretical significance), but a useful tool as well in the everyday 
practice of working relativists. The successes in the proof of the positivity of energy inspired several relativists to 
search for expressions of the gravitational energy-momentum at the quasilocal level, i.e. to associate these physical 
quantities with closed spacelike 2-surfaces 5 ( see, for instance, Ref. [5] and the references therein). 

Quasilocal energy is important for its several potential applications. It can be used to define binding energy of 
stars. In numerical relativity it can tell us how to cut off the data on a noncompact region to a compact region where 
the quasilocal mass of the compact region approximates the ADM mass of the noncompact region. When we evolve 
the data according to Einstein equations, we need local control of energy to see how the space changes. The positivity 
of quasilocal energy is essential for such an investigation. This could be regarded as a generalization of the energy 
method in the theory of a nonlinear hyperbolic system. 

There are several contenders for a good definition of quasilocal energy. The two that interest us are Brown and 
York’s ’canonical quasilocal energy’[6], and the various definitions based on the integral over S of the Witten-Nester 
two-form (the two-form used in Witten’s proof of the positive energy theorem [3]). The approach of Brown and York 
starts with standard Hamilton-Jacobi theory and thus is somewhat different from more traditional approaches based 
on Noether techniques. Brown and York have determined what geometric entity play the role of quasilocal energy in 
general relativity. The canonical quasilocal energy is supposed to be the energy of the gravitational and matter fields 
contained in a spatial volume E, whose boundary two-surface is S = <9E. This approach gives the correct quasilocal 
energy surface density and allows one the freedom to assign the reference point of the quasilocal energy. 

The approach based on the Witten-Nester two-form uses spinorial methods, and the different definitions are distin¬ 
guished by the choice of supplementary equation the 5-spinors are supposed to satisfy, for example the Sen-Witten 
equation [3], the Dougan-Mason equation [7], or the Ludvigsen-Vickers equation [8]. Since Witten’s positive energy 
proof can be understood in terms of the Hamiltonian, the Hamiltonian density associated with this proof provides a 
locally positive localization—and thus has real promise as a truly physical energy-momentum density for gravitational 
fields. To fulfill this promise certain features need further consideration; an outstanding one concerns the role of—and 
even the need for—the spinor field, which seemed rather mysterious. Although there are very beautiful argument for 
some nice important results, it is not clear as to how much the Witten argument really captures the correct physics. 
Moreover, all the spinor formulations cannot give angular momentum and the center-of-mass moment [9]. 

The two approaches are conceptually distinct, even though in some circumstances one might expect their numerical 
value to coincide. Therefore, it is of interest to establish a connection between the Brown-York quasilocal energy and 
the Witten-Nester form. Indeed, such a connection supports the use of the Witten-Nerster form in the quasilocal 
context, and, furthermore, provides one with a non-spinor vantage point from which the various supplementary spinor 
equations can be examined. The general relationship between the Brown-York quasilocal energy and the common 
’’spinorial definition of quasilocal energy” constructed from the Witten-Nerster two-form has been established by Lau 
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[10]. He found that the Sen-Witten energy expression can always be written in a form which is essentially identical to 
the that of the Brown-York quasilocal energy. His theory is an Ashtekar-variable reformulation of the metric theory 
of quasilocal stress-energy-momentum originally due to Brown and York. It is meaningful to apply Lau’s formalism 
to the proof of the positive energy theorem. In this paper we start with the self-dual Lagrangian following Lau’s route 
and then use Witten-Ashtekar-Horowitz’s method [11] to prove the positive energy theorem. To this end we have to 
derive a set of extended Witten identities and then apply them to the positivity proof. In the discussion on positivity 
proof and energy-momentum quasilocalization the boundary term has important roles. The value of the Hamiltonian 
”on shell” is determined by the Hamiltonian boundary term. The ’’natural” Hamiltonian boundary term inherited 
from the Lagrangian can—and should be adjusted. The general expression for boundary term depends on the choice 
of variables, a displacement vector field (e.g. translation for energy momentum and rotation for angular momentum), 
a reference configuration and boundary conditions. Being different from most of the works on positivity of energy, 
in this paper by choosing a boundary term with a nonvanishing shift we will extend the proof to involve momentum 
and even angular momentum. 

The organization of this paper is as follows. In Sec. II. a review of some results of Lau is given in a rather different 
appearance. The Lagrangian of general relativity is written in terms of two-spinors and then is decomposed into 
two independent parts, the self-dual and the anti-self-dual part. The self-dual part is identified with a version of the 
Ashtekar Lagrangian [12] given by Goldberg [13] and then can be taken as the gravitational Lagrangian. Sec. Ill is 
the Hamiltonian formulation in which the dyad spinors are chosen as the basic dynamic variables while the Ashtekar 
connections appear in the conjugate momenta. In Sec. IV the Hamiltonian boundary term is rewritten as a suitable 
form for the application of the extended Witten identities. The extended Witten identities are derived in Sec. V and 
are used to calculate the boundary term and prove the positivity of the energy in Sec. VI. Finally, Sec. VII is devoted 
to some conclusions and remarks. 


II. THE SELF-DUAL LAGRANGIAN OF GR 

We start by rewriting some results of Lau [10] in a rather different appearance in order to adapt to our purpose. 
Using the spin connection one-form u> a KL the curvature two-form R a bKL of the spacetime M can be written as [14] 

RabKL = ^a^bKL — V b ~ U a M K^bML + K^aML, (1) 

where is the covariant derivative defined by the affine connection T a bc , for example, the covariant derivative of 
the cotetrad one-form e J b is 

V^e J b = e\u al J = d a e Jb - T c ab e J c . 

Then the scalar curvature R can be written as 

p _ 0X7^ (£>^ a \ I Ka Lb I Ka Lb I /o\ 

K — zv a e LObKL) + e e cu a k^ml ~ e e u a L WbKi, (*) 

and the gravitational Lagrangian Lq consists of the Moller Lagrangian [15] and a divergence term: 

£g = V-gR 

= V~9 (e Ka e Lb u) a I K^biL - e Ka e Lb UaL I u b Ki) + 2<9 a ( V zz ge Ka e Lb uj bK L) ■ 

Passing on to two-spinor expression [16], the spin connection reads 

CC' C C' _ Cf c 

U AA'BB' ' = WAA'B C_B' + UAA'B' 

_ , CC' ,,,- CC' (o\ 

^ AA'BB' TW AA'BB' ’ 

where 

, ,+ CC' _ , , c c c' ( A \ 

10 AA'BB' - UAA'B €B' j. (4) 

and 

— CC' — Cc C (tr\ 

A> AA'BB' — UAA'B’ £B , (p) 

is the self-dual and the anti-self-dual part of the connection u>aa'bb' CC , respectively, and then the scalar curvature 
R splits into two parts 
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R = R + + R~, 

where the self-dual and the anti-self-dual part of R is, respectively 

p+ OY7^ ( C AT K'a LZ/ b \ i K I ,LK' KK' I L 

R = 2V a le e UJbKLtL'K' I + U> K> /fW 7L - OJ LU K'lK 

d- _ nvy r („KK'a 0 LL'b— , \ , — A"/' — LL' —KK'I' — L> 

K — zv a I e e cobK'L'eLK J + wj? at'W bl'—w l'^k bk’- 
Consequently, the gravitational Lagrangian Cq = yJ—gR splits into two pats 

Cg = C + g + C~ g . 

The self-dual part Cq consists of the self-dual Moller Lagrangian [13,15] and a divergence term: 

r+ ( 4 ) ^ (, ,K I ,LK' , KK'I , L \ (n(4) KK ' a LL 'b , A 

t-G ° l W K> K<jJ IL—W LU K'lK) + da ere e UbKLtL'K'l 

where 

(4) a = V=? 


( 6 ) 

( 7 ) 

( 8 ) 


(9) 


( 10 ) 


is the determinant of the SL(2,C) soldering form cr fl AA on the spacetime manifold M. Since Cq and C G depend 
on u>aa'b C and ZJaa'b ,Ci respectively, and are independent of each other, we can choose the self-dual part Cq as 
the Lagrangian of the gravitational held, which is identified with a version of the Ashtekar Lagrangian [12] given by 
Goldberg [13]. 


III. THE HAMILTONIAN FORMALISM OF THE SELF-DUAL GRAVITY 

To put the theory in Hamiltonian form, we will assume that the spacetime M is topologically E x R for some 
spacelike submanifold E and assume that there exists a time function t with nowhere vanishing gradient {dt) a such 
that each t =const. surface E t is diffeomorphic to E. t a will denote the time flow vector held ( t a ( dt) a = 1), while n a 
will denote the future-pointing timelike vector held ( n a n a = —1) normal to the t =const. surfaces. Then the intrinsic 
metric q a b of E, the lapse N , and the shift N a are induced by 


9ab Qab ei a Tlb, 


and 


t a = Nn a + N a , 

respectively. In the spinor notation [16], the normal vector n AA dehnes an isomorphism from the space of primed 
spinors to the space of unprimed spinors: 


x+ A = iV2n AA ' 


OLA' 


UABCD = 


iV2n B A A) A A' 


CD , 


In this notation 


n 


AB = is/2r BA '- A 


n a' 


i AB 

V2 


( 11 ) 


Qab — Qab '^a'^bi 


reads 


gABCD = q ABCD ^ n AB n CD^ 


or 
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( 12 ) 


e 


AC e BD 


_ e A(C e D)B 


}_ AB CD 
2 


In terms of the dyad spinors 


(o C = o c ,(i C = L C , 


the self-dual connection can be expressed as 


vabc D = Cc b V r AB t b D , (13) 

Using these results we can decompose the quantities appearing in Cq into the tangential and the normal parts to the 
surface E t , for example, for a four-vector V AB we have 

V AB = q AB V CD - n AB n CD V CD 
= y( AB ) y\ AB \^ 

where V (AB) and V [AB] is the tangential and the normal part of V AB to the surface £ t , separately. Then the 
Lagrangian Cq can be rewritten as 

£+ = As/2icj^ AB \ a B C aA 

-Na[V r (AB) C A a V r ^ CB \c a - v[ AB) Cc a V r ( sc )C A a] 

-y/2iaN CD (3V^ D a a V r(AB) CaB + (A a V£ D V r(AB) (as) 
d (AB) (2Na( Ab V r ( CB )( bC - V2iaN CD ( Aa V r {CD) (a B ) , (14) 


where 


and 


CaA~ t CD Vc D (aA, 


= 4 _z = VEl 

N N ' 


(15) 


(16) 


The dyad spinors £ a A are chosen as the basic dynamic variables and then the canonical momenta conjugate to ( a A 
are 


~aA = d£G_ = 4 ^ 2 iaV T(AB) c 


d(aA 

The gravitational Hamiltonian density can be computed 


n G =p aJl ( aA -c G 

= ^[v[ AB) C A aV r ( CB >Cc a - Vf AB) Cc a V r ( BC )C A a] 

+ V2iaN CD [3 (V£ D Ct°) V r ^ S )CaB + (A b V r CD V r(AB) ( b B 
+d {AB ) (2NaC Ab ^ CB \ bC - V2iaN CD ^ Aa W T (CD) Ca B ) . 


(17) 


(18) 


( 19 ) 


Using 




~caA 


4V2 0 


we obtain 
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^ = ~ aV UB)Cc a ^ r(BC \ A a} 


--N CD ^3 (VZnCaA) V aA + < aA b S7 T CD {-F 
+d { AB) [^NQ a A F B - V2iaN CD ( Aa V r (C D)Ca 


( 20 ) 


i. e. 


where the Hamiltonian constraint 


H g = NH + N ab Hab + d {AB) B 


the momentum constraint 


U = 35a PaAl?aA ~ aV UB)(*cV r(BC) C A , 


( 21 ) 

( 22 ) 


and the boundary term 


Hab = -7 F^bU - ^(aC^AB) ( ) , 


B ab = 2Na^ Ab S/ T{CB) C, bC - V2iaN CD C Aa V r (CD) Ca B 


s/2i 


NCa A p aB - V2iaN CD C Aa \7 r (CD) C 


B 

a ? 


(23) 


(24) 


separately. Being different from Ashtekar’s formulation, here the dyad spinors Q a A are chosen as the basic dynamic 
variables and the conjugate momenta p aA are related to Ashtekar’s variables oj^ab) CD [12,13] by 

“ {ac)bC = IT^ Cb “ ? “- 4 ' 


IV. BOUNDARY TERMS 

Being different from most of the works on positivity of energy in which the shift N AB is chosen to vanish, we 
consider a more general boundary term B ab which allows a particular choice of N AB . Choosing the triad 


_ AB * (AB ,A,B ' : 

a i = —= [oo —Lb 

V 2 


) = ~j= (™ a + m a ) 


*2 AB = ^=(o A o B + L A L B )=-^=(m»- 


V2 

a3 AB = j_ ( _ 0 A l B _ l A o B 
v 2 

on £ then the normal unit vector to the boundary <9E is 


V2 


m 




(25) 


V““ = <J 3 “" =- -= ( o a l b + l a o b ) . 

V2 v ' 


We chose 


N = -y 2 , 

N AB = X 2 (m AB + m AB +v ab ) 

= -^=x 2 (o A o B + l a l b - o^ a l b) 


(26) 

(27) 
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and compute the boundary integral 



dSAB 


B AB v AB dS 


— cry 2 (m b m a Vb/ a — m b m a Vbn°) dS 

£ a X 2 (j n>> ~ mb - p b + \ n>> ^j ( n a^bl a - m 0 V & ?n°) dS 

- j) a\ 2 k - ^a\ 2 (m b - m b - | l b + (n 0 V^° - m 0 V b ?n°) dS 


(28) 


(29) 


where k is the trace of the extrinsic curvature of <9£. The first term of the integral is just the ”unrefernced” Brown- 
York quasilocal energy when x 2 — 1 [6,10], the second term comes from the particular choice (27) of N AB and then 
includes the contribution of momentum and angular momentum on which our interest is concentrated. 

For further calculation and discussion we rewritten the expression of B ab as follows. 

Since 

d(AB) ('V2iaN CD ujcD AB ) 

= V2iaDf2 B) ( N cd lo CD AB ) 

= V2iae AE e BF e CG e DH D^2 f) (Ngh^cdab) 

= V2ia ( e DB e FH + e FD e BH ) (e CA e EG + e EG e AG ) D$ F) ( N gh u C dab) 

= V2iaDf2 B) (- N ab u cd CD ) - V2iaDf2 B) (- N ad u c BC d) 

-V2iaDf2 B) ( N ca u bd CD ) + V2iaDf2 B) ( N cd u ab C d ) 

= V2id {AB ) (aN AB ( Cb V r {CD) ( b D + aN GD Cc b ^ rAB Cb D ) , 

B ab can be written as 

B ab = 2Na^ Aa V T{CB \ aC - V2iaN AB C Cb V r {CD) C b D 



~V2iaN GD Cc a ^ nAB \aD. 

(30) 

Introducing spinors A A 

and A +yi by 



A 1 -v A A 1 \+A 

o = —A , i = —A 

(31) 


X X 

and noting 

V r AA ’\ C = d A A>\ C i 

(32) 

we can compute 



B ab = 

2NaX (a a 9 (cb) AJ - A+^)Ac) + 27Vo' —(A +a A c - A a AJ) 9 (cb) x 
-V2iaN AB X (X C d {CD )X +D - X +c d {CD )X D ) 

X 

-V2iaN CD X (Ac5^A+ - X+d^X D ) . 


Noting 

N = -x 2 = -Xa\ +A . 

(33) 

and 
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(34) 


N ab = -m AB - m AB - v AB 

= -J= (x A x B + X +A X +B + A (a A + b) ) , 

we obtain 

B ab = - J (a a A b + A +a A +b + A (A A +b) ) {X c d {CD )^ +D - A +c d (CB) A B ) 

-2a (x A d^X + c - X +A d^xc) - a {X +A X C - X A X + C ) 8^ ln X 2 
+a (V B d (AB) A+ + X D d {AB) X D + ^X D d (AB) X+ + ^X +D d {AB) X d) ■ (35) 


V. THE EXTENDED WITTEN IDENTITY 


In order to deal with the boundary term B ab we follow Witten [3], Ashtekar and Horowitz [11] and derive some 
useful identities. In a two-spinor formulation a Dirac spinor e and Dirac matrices j 1 can be expressed 


, A 

£ = I - 

TA' 


/ = V2 


0 a 

a 1 a'A 0 


IAA' 


respectively, and then the Witten identity [3] 

i 

-V 2 £ = (iV) 2 e = yi v - + 4 ^ G ( T oo + 5] T 0j 7 y>. 

i j 

becomes two two-spinor identities 

-2V ag V GB A b = -V bg V bc A a + 47rG(T 00 A A + V2T 0 A B X B ), 
— 2Va'C'V t b ’ = -V 


(36) 


(37) 


'X7b'C' t a' + 4ttG(Toota’ + V2Tqa ,b tb>), 


where V ab and V a'B' denote the covariant derivatives defined by the spin connections ojabc D and u>a'B'C' D 
respectively. Multiplying the first identity with A +a leads to 


Since 


2A +a V (AG) V (CB) A b = A +a V (BC) V( bc )Aa - 47tG(TooA +a Aa + V2T 0 abX +a X b ). 

a +a v (BC) v( bc >a a = v (BC) (a +a v( bg >a a ) - (V (bc) A + a ) v^ BC h A 
= D bc (A +A V( BC )A a ) - (V (bc) A a ) + V (bc U^, 

then we have 

aD B c (a +a V( bc >A A ) = 2aX +A \7 { AC)V iCB) X B + a (V (BC) A a ) + A^ 

+4^Ga(T 00 A +A A^ + V2T 0AB X +A X B ). 

Integrating it leads to 

j> <jX +A ^ BC h A dS { AB) = j <tD bc (a +a V ( bg >A a) dv 

= 2 f aX +A V { AC)V {CB) X B dV+ [ a(y {B C)^ A ) + ^ {BC) X A dV 
J s J s 

+47tG J <7(T 00 X +a X a + V2T 0AB X +A X B )dV. 


(38) 
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Multiplying the first identity of (37) with A' 4 leads to 

Dbg (a a V (BG) Aa) = 2A a V ( ac)V (GB) A S + 4,V2ttGT 0AB X A X B . 

Integrating yields 

j <j\ A V (BC) \ A dS (AB) = j oD bc (a a V (bg Ua) dV 

= 2 [ <r\ A \7 (AC) V {CB) \ B dV+ 4V2irG f aT 0AB X A X B dV. 

JY J Y 

By the similar way, from the second identity of (37) we obtain 

2 [ a\ A 'v {A , cl) ^ C ' B 'h B/ dV + AV2nG f a\ A ' T 0A , B '\ B ,dV 
J Y JY 

-n/2 j aX +B K CD A B X A dS iBC) 


aX +A V a XtdS a = 


and 


j crX A \7 a X\dS a 


A '^(A'C')^ C ' B,) X B 'dV + J V (BC) A a 


= 2 g A Vm'cmV 
J Y 


(v( bc >a 


■a I dU 


+47tG J ct(T 00 X a ' X a > + V2X A 'T 0A , B 'X B ,)dV + iV2 j> aX B K CD A B X\dS {CD) 
-iV 2 J oK C da B (a a V (GB >A+ - A+V (GB U A ) dV + 2j gK C da 


E K CDB E X A X+dV. 


(38), (39),(40), and (41) can be called the extended Witten identities. 


(39) 


(40) 


(41) 


VI. THE CALCULATION OF THE BOUNDARY INTEGRAL 


Now we are ready to calculate the boundary integral. Since 

1 . 1 

X ' X 


CoA — OA — —X A, Cl A — L4 — — A + 


A> 


then 


and 


V^,A g = d AA 'X c , 


UAA'B — Cb V AA'Cb 


= 4 (Ab£WA +c - A+cWA c ) + 4^AA'X (A£A g - A b A + g ) . 


X 


Vaa'Ab : = d AA 'X B — u AA < B Xc 


= d A A'X B - 2 (Ab9aa'A +c - Xgd AA >X c ) Xc H— Oaa'XXb- 


X 


X 


Using this result we obtain 


A +s Vaa'A b = 2X +B d AA 'X B - xQaa'X, 


and 



\ b V AA '\b = X B d AA> X B - d AA 'X c Xc = 2X B d AA >X B . 


By the same way we find that 


which leads to 


Then we have 


X7aa' Ajj — d AA i Ag u) AA 'b A^ 


= d AA ,X+ - (X B d AA ,X +c - X B d AA >X c ) X+ + -d AA , X Xtn 


A VaA'^b ~ ^A QaA'^b + X^AA'Xi 


X +B V AA ,X+ = 2X +B d AA> X+. 


Oaa’^b H~ Daa'^b — 2 (A +5 Vaa'A_b + A^Va^'A^) 
X B d AA >X B = -A S V aa >X b , 

X +B d AA ,X+ = ^X +B V AA ,X+, 


X+ D d ( AB )x D + X D d^X+ = i (a +b V^Ab + A 1 


x D d^x D = -a b v( ab Ub, 

2 

a +o (9 (ab)a+ = I A +b V (j4b) A+. 


Using these results and (35) we obtain 


j> B AB dS {AB) = ~f^2 (X A X B + X +A X +B + A^A +S )) (A c 9 (cd) A+ b - A +G a (C B)A B ) dS {AB) 

—2<j><7 ^X A d^ CB ^X^ — A +yl <9 (CB ^Ac^ dS( AB ) 

- j> ° ^ (A +a A c - A a AJ) d (CB) x + |<9 (AB) x dS {AB) 

T~ (j) & (V b V<9 (ab) A+ + A d V (ab) A d + a +b V (ab) A b ) dS^s) 


From the extended Witten identities (38), (39), and (40) one finds 


i f a (V B V (AB) A+ + A B V (AB) A B + A +B V (AB) A D ) d5 (iB ) 

[ <jX +A \7 (AC) V (CB) X B dV + [ aX A \/ (AC) V {CB) X B dV+ [ ctA A 'v (a , c 0 V (C ' B ^A B -dU 
Je is is 

+ lt a ( V (bc)A A ) + S7 {BC) X A dV + 2nG j^a(T 00 X +A X A + V2T 0AB X +A X B )dV 


aX +B K CDA B XtdS l 


B A A ad(CD), 


and then the boundary integral becomes 
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f B AB dS {BC) 

- 4 (a a a s + A +a A +b + A< a A + b )) (X c d {CD) X +D - A +c a (CD) A D ) ds (AB) 

—2 f a (Vd (GB) Aj - X +A d { - CB) X c) dS(BC) 

~ f a x ( a+Aa c - A a A£) <9 (gb) x + ^x<9 (AB) x c!S(bc) 

+ [ aX +A \7 (AC) V {CB) X B dV + f aX A V {AC) V^ CB) X B dV+ f crA^V^qV^'^A B 'dV 
J s ./E JE 

+ lL a (y(BC)X A ) + V {BC) X A dV + 2nG fa(T 0Q X +A X A + V2T 0AB X +A X B )dV 
~~Y' / aX+BKCDA B^A d S(CD)- 


When the Witten equations 


V (GB) A B = 0, V^ C ' B '^X B ' = 0, 


are satisfied we have 


Noting 


f B AB dS (AB) 

- 4 ( aAaB + a+Aa+B + A( a A + b )) (A g 9 (cd) A+ b - A +G d (CD) A B ) ^ (AB) 

J sC / 

—2 f a (^X A d < ' CB ' ) X q — X +A d l ' CB )XdS( AB ) 

-fa l(X+ A X c -X A X+)dW X +\xdt AB h dS (AB) 

+ \j a {^(BC)X A ) + ^ {BG) X A dV + 2nG fa(T 00 X +A X A + V2T 0AB X +A X B )dV 
/ aX+BK ° DA BX\dS(CD)- 


dS( AB ) = V( AB )dS, 


one finds 


AB 1 ( A B , A B\ 

V = - -= [O L + b O J , 

v 2 

/ J ( aAaB + a+Aa+B + a(Aa+B) ) ( AC a (CD) A +D - X +c d {CD )X D ) dS (AB) 

—= df cr (A g < 9( C - £) )A +j:> — A +g 3(C£))A d ) dS. 


And according to the definition of the extrinsic curvature of £ 

KcDAB = —V(CD) a 0AB 


= -fjf(CD) {-0 A L B + L A 0 B ) 


we have 
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X +B KcdabX +a = 0. 


Finally, the boundary integral becomes 


B ab dS(AB) = ° (^ C ®(cd)X +d — X +c 8(cd) X D ) dS 

—2 y er ^A' 1 2 3 4 5 6 7 8 9 10 11 9 < ' CB ^A[t — A +a 9^ c ' b ^Ac^ V(AB)dS 


/ ^ [ A ^ A+AAc ~~ A ^) d{CB)x + \ xd(AB) X v (AB)dS 

l [ <j(W (B c)X A ) + X/ (BC) X A dV 


+2t tG J a(T 00 X +A X A + V2T 0 A B X +A X B )dV. 


(44) 


This is manifestly non-negative if the gravitational source satisfies the dominant energy condition and then the 
positivity of the quasilocal energy can be proved easily. Starting from the choice (34) of the shift vector N AB , the 
prove can be extended to include momentum and angular momentum. 


VII. CONCLUSIONS AND REMARKS 

Using a spinor method the Einstein-Hilbert Lagrangian is decomposed into the self-dual and the anti-self-dual part. 
The self-dual part is identified with a version of the Ashtekar Lagrangian given by Goldberg. Starting from this 
Lagrangian and choosing the dyad spinors as the basic dynamical variables a Hamiltonian formulation is developed. 
A set of extended Witten identities is derived. Choosing a particular shift vector N AB and applying the extended 
Witten identities the proof of the positive energy theorem is extended to a case including momentum and angular 
momentum. 

In both the approach of Brown-York and the approach of Witten-Nester the boundary term has important roles. 
The ’’natural” Hamiltonian boundary term inherited from the Lagrangian can—and should be adjusted. The choice 
of the Hamiltonian boundary term is the key point in the discussion of this paper. Being different from most of 
the works on positivity of energy in which the shift N AB is chosen to vanish, we consider a more general boundary 
term B ab which includes the ’unreferenced’ Brown-York quasilocal energy and allows a particular choice of N AB . 
However, we do not deal with the issue of reference point since it has been shown [10] that when the Witten-Nester 
expression is evaluated on solution spinors to the Sen-Witten equation (obeying appropriate boundary conditions), 
an implicit reference point for the energy is set. The further investigation of the boundary term in some particular 
cases will be given in a forthcoming paper elsewhere. 
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